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An attempt is made to clarify the difference between a the- 
orem derived by Evans and Searles in 1994 on the statistics of 
trajectories in phase space and a theorem proved by the au- 
thors in 1995 on the statistics of fluctuations on phase space 
trajectory segments in a nonequilibrium stationary state. 
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Recently a Fluctuation Theorem (FT) has been proved 
by the authors (GC), [Q, for fluctuations in nonequilib- 
rium stationary states. Considerable confusion has been 
generated about the connection of this theorem and an 
earlier one by Evans and Searles (ES), 0, so that it 
seemed worthwhile to try to clarify the present situation 
with regards to these two theorems. 

In a paper in 1993 by Evans, Cohen and Morriss ||, 
theoretical considerations lead them to a computer exper- 
iment about the statistical properties of the fluctuations 
of a shear stress model (viscous current) - or the related 
entropy production rate - in a thermostatted sheared vis- 
cous fluid (plane Couette flow) in a nonequilibrium sta- 
tionary state. 

The Fluctuation Relation found in the simulation 
reads in current notation, [flj: 
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Here 7r T (p) is the probability of observing an average 
phase space contraction rate (which in the models con- 
sidered has the interpretation of average entropy produc- 
tion rate) of size p a + on one of many segments of dura- 
tion r on a long phase space trajectory of the dynamical 
system modeling the shearing fluid in a nonequilibrium 
stationary state; here cr(x) will denote the phase space 
contraction rate near a phase point x (i.e. the divergence 
of the equations of motion) and <r+ is the average phase 
space contraction rate over positive infinite times so that 
p is a dimensionless characterization of the phase space 
contraction (with time average 1). The approximation 
within which eq. ([!]) was observed was very convincing 

§• 

Under suitable assumptions, see below, a more precise 
formulation of eq. ([!]) was derived in Q| : 
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where the validity of the fluctuation relation eq. ([!]) for 
asymptotically long times r is more clearly expressed. 
Later several other computer experiments have confirmed 
the relation eq. (§),!,[§,§._ 

The original computer experiment, g], was inspired 
by a theoretical argument for the relative probabilities 
to find a phase space trajectory segment of length r in a 
state x with phase space contraction rate p and in a state 
x' with rate —p. These theoretical considerations lead to 
the correct prediction eq. ([!]), which was confirmed by 
the - independently of the theory - carried out computer 
experiment. 

Evans and Searles gave a derivation of a theorem 
which had a similar form as eq. (|l]). More precisely: let 
E p be the set of initial conditions of a dynamical system 
for phase space trajectories along which the phase space 
contraction is e~ pCT+T in a time T. We denote by /j,l(E p ) 
its Liouville measure. Similarly, let (1l(E- p ) be the Li- 
ouville measure of the corresponding set of phase space 
trajectories along which the phase space contraction in 
time T is e p<T+T . Quite generally, and in all models con- 
sidered in the literature relevant here, E_ p = IStE p , if 
St is the time evolution (Liouville) operator of the sys- 
tem, and if I denotes the time reversal operation, so that 
t — > StX is the phase space trajectory over time t start- 
ing at x at t = 0. Hence E- p , the set of points around 
which phase space contracts at rate — pcr+T ', is obtained 
by evolving forward over a time T those in E p (which 
would contract by pa+T) and then inverting the veloci- 
ties by the time reversal operator I. In fact, the sets E p 
and E- p are those considered by Evans and Searles in 



Then the proof in j| is the following: 
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where one has used that the Liouville distribution is time 
reversal invariant, i.e. pr,(E) = fiL(IE) (although it is 
not stationary) to get the second equality as well as the 
definition of phase space contraction in the third equality. 

The arbitrary time interval T includes the short times 
referring to the transient behavior of the system be- 
fore reaching the nonequilibrium stationary state. In 
the derivation of eq. (JsJ) only time reversal symmetry 
is used. Later, [Q, it was argued that under this as- 
sumption alone, eq. (^) also holds in the nonequilibrium 
stationary state /Xoo> since eq. (||) is valid for any T and 



1 



the Liouville distribution /i^ would evolve in a sufficiently 
long time T into a distribution St^l arbitrarily close to 
a nonequilibrium stationary state fi^. 

Therefore the eq. (|J) was reinterpreted in and as- 
serted to be identical to eq. (|^), which, however, refers to 
the statistics of trajectory segments, along a trajectory 
in a chaotic nonequilibrium stationary state /ioo, not to 
the statistics of independent trajectory histories emanat- 
ing from the initial Liouville distribution fiL under the 
time reversibility assumption. 

In 1995 the authors proved eq. (j^) based on a dynami- 
cal assumption, called Chaotic Hypothesis (CH), which 
assured strong chaoticity ("Anosov system-like behav- 
ior") for the systems for which eq. (Q) held. In that 
work the name Fluctuation Theorem (FT) was first in- 
troduced for eq. @ , and was proposed as an explanation 
for the experimental result eq. (|l|). We will call this the 
GCFT. 

It is worthwhile to emphasize again that, while the 
right hand sides of the eqs. (0) and (§) look very similar, 
they, as well as the left hand sides of these equations, 
really refer to entirely different physical situations. 

The eq. (^), 0], holds for any T on trajectories with 
initial data sampled from the Liouville distribution at t = 
and it can be considered as a simple, but interesting, 
consequence, for reversible systems, of the very definition 
of phase space contraction. We will call it here the ESI, 
where the I refers to "identity". The ESI is much more 
general than the FT in eq. (^), which needs, in addition 
to phase space contraction (<r + > 0) and time reversal 
symmetry, also the Chaotic Hypothesis. The proof of 
the ESI, fully described in eq. (|3|) above, is identical in 
essence to the proof in Q which is much more involved. 

In order to illustrate the fundamental difference be- 
tween the two theorems we first give an example of a 
case where the more general ESI eq. @ holds, while the 
GCFT eq. (|) does not. 

To that end we consider a single charged particle in 
a periodic box, with charge e moving in an electric field 
E, i.e. a Lorentz gas without scatterers, and subject 
to a Gaussian thermostat (to obtain a nonequilibrium 
stationary state): 

q = p, p = eE ap (4) 

where the "thermostat" force — ap, with a — as- 
sures the reaching of a nonequilibrium stationary state 
of this system. 

In this case one can solve explicitly the trivial equa- 
tions of motion eq. (|J) and check that the Liouville dis- 
tribution (1l indeed evolves towards a stationary state 
/ioo, which is simply a state in which the particle moves 
with constant speed parallel to E. The ESI eq. (||) will 
(of course) hold for the phase space trajectories of this 
system sampled with the initial Liouville distribution /i£, 
but it will not be a fluctuation theorem, since there are 
no fluctuations. Also, GCFT's eq. (||) will not hold for 
the phase space trajectory segment fluctuations of this 



system, which is not a contradiction because the system 
is not chaotic. 

From this simple example follows that the two theo- 
rems cannot be equivalent, and the validity of eq. (Q) 
cannot imply much, without extra assumptions, on the 
fluctuations (absent in this case) in the stationary state. 
Note that eq. (|) is an identity which is always valid in 
the systems considered. The system of eq. (|j) is there- 
fore a counterexample to the statement, that eq. (||) im- 
plies eq. (||), i.e. to the statement, (?]], that ESI implies 
GCFT. 

Second, and more interestingly, one can try to derive 
the GCFT eq. (Q) from the more general ESI eq. ||). 
One could then try to proceed as follows. 

First one would need to show that on a subsequent 
trajectory segment of length r, after time T, the ratio 
of the probabilities of finding a phase space contraction 
of +pa + r to that of finding — per+r over this segment of 
length r, would be given by e pa+T . Here po+r is any 
preassigned value of the phase space contraction. How- 
ever, eq. ^ gives no information whatsoever about those 
points which after a time T evolve into points which in 
the next r units of time show a phase space contraction 
±ptr+T. In other words, the ESI does not contain the 
detailed information needed to derive the GCFT. 

If one adds the Chaotic Hypothesis to the time re- 
versal symmetry assumptions made about the dynam- 
ical system in the ESI, one could use Sinai's theorem, 
, to assert that such a system, starting from the ini- 
tial Liouville distribution /i£, will indeed approach a 
chaotic nonequilibrium stationary distribution /i^ sup- 
ported (however) on a fractal attractor A with Liou- 
ville measure hl{A) — 0. This is, in this case, the SRB 
distribution, fisRB , of the system, which was used in [[l] . 
However, for a proof of the GCFT, details of the SRB 
distribution are needed, which contain just the details 
considered in That is, one has to make an appro- 
priate (Markov) partition of the phase space and assign 
weights to the cells of increasingly finer partitions lead- 
ing, to the SRB distribution. This then allows one to 
assign appropriate weights to those regions (of Liou- 
ville measure) in phase space that will give rise, to phase 
space contractions on trajectory segments r of ±po~ + T 
leading to the GCFT. 

A more advanced comparison between the ESI and 
GCFT requires a more quantitative statement of the lat- 
ter result. Namely, 0], eq. (||) can be derived from the 
stronger relation: 

MP) e ( P . + + o(^))T (5) 

where Too is a time scale of the order of magnitude of 
the time necessary in order that the distribution St^Li 
into which the Liouville distribution (i^ evolves in time 
T, be "practically" indistinguishable from the stationary 
state that we denote by /ioo. Here the validity of the 
fluctuation relation (Q) for asymptotically long times t 
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is more clearly expressed. The esistence of the time 
and its role in bounding the error term in eq. (^|) are 
among the main results of fiL The time appears 
in jj] as the range of the potential that generates the 
representation of 

Moo cl Gibbs state using a symbolic 
dynamic representation of the SRB distribution on the 
Markov partition, || . 

Examining the ESI derivation above, one easily sees 
that the following 

StHl(E p ) = fi L (S- T E p ) = c ( pj++0 (Z)) T ^ 

can also be derived, p| ; it is important to note that 
the argument leading to eq. (|^) cannot say more than 
this: in particular one must justify why T, which in prin- 
ciple should be large, strictly speaking infinite, so that 
St^l be identifiable with fJ-sRBi can in fact be taken 
smaller than r. Justifying this requires assumptions (as 
the above counterexample indicates) like the mentioned 
Gibbs property of the SRB distribution which is even 
stronger than requiring exponential decay of correlations 
(also implied by the CH). Otherwise, without extra as- 
sumptions, one has to say that T has to be taken infinite 
first with the result that eq. (||), hence eq. (||), becomes 
empty in content. Of course one can argue that "on phys- 
ical grounds" T needs not to be taken infinite but just as 
large as some characteristic time scale for the approach 
to the attractor: but the precise meaning of this, and 
the assumptions under which it can be stated, is pre- 
cisely what needs to be determined particularly because 
in nonequilibrium systems the attractor A is fractal with 
/ii(A) = and one can very well doubt that the StI^l 
distribution is ever close enough to the /loo distribution 
to allow comparing eq. (^|) and (^). This requires a con- 
vincing argument, were it only because St^l{A) = 0. 
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